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Abstract 

The purpose of this paper is to present simple and general alge- 
braic methods for describing series connections in quantum networks. 
These methods build on and generalize existing methods for series (or 
cascade) connections by allowing for more general interfaces, and by 
introducing an efficient algebraic tool, the series product. We also 
introduce another product, which we call the concatenation product, 
that is useful for assembling and representing systems without neces- 
sarily having connections. We show how the concatenation and series 
products can be used to describe feedforward and feedback networks. 
A selection of examples from the quantum control literature are ana- 
lyzed to illustrate the utility of our network modeling methodology. 

Keywords: Quantum control, quantum networks, series, cascade, 
feedforward, feedback, quantum noise. 

1 Introduction 

Engineers routinely use a wide range of methods and tools to help them 
analyze and design control systems. For instance, control engineers often 
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use block diagrams to represent feedforward and feedback systems, Figure 
[TJ Among the methods that have been developed to assist engineers are 
those concerning the connection of components or subsystems to form a net- 
work. One of the most basic connections is the series connection, where the 
output of one component is fed into the input of another, as in Figure [TJ 
When the components are (classical, or non-quantum) linear systems, the 
connected system can be described in the frequency domain by a transfer 
function Q{s) = G2( s )Gi(s) which is the product of the transfer functions 
of the components. The description can also be expressed in the time do- 
main in terms of the state space parameters Q = (A, B, C, D) (as we briefly 
review in section [2]). The series connection has an algebraic character, and 
can be regarded as a product, Q = Q 2 < G\- Because of new imperatives 
concerning quantum network analysis and design, in particular, quantum 
feedback control, [21], [25], [H], [23], [26], g], pi], [12] the purpose of this 
paper is to present simple and general algebraic methods for describing series 
connections in quantum networks. 
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Figure 1: Series connection of two (classical, or non-quantum) linear systems, 
denoted Q = Gi < Gi- 

The types of quantum networks we consider include those arising in quan- 
tum optics, such as the optical network shown in Figure [21 This network 
consists of a pair of optical cavities (discussed in subsection 13. 2j) connected 
in series by a light beam which serves as an optical interconnect or quantum 
"wire" . In this paper (section we show how series connections of quantum 
components such as this may be described as a series product G = G 2 < Gi. 
This product is defined in terms of system parameters G = (S, L, H), where 
H specifies the internal energy of the system, and X = (S, L) specifies the 
interface of the system to external field channels (as explained in subsection 
13.41 and section |4j . 
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Figure 2: Series connection of two optical cavities via an optical interconnect 
(light beam) or quantum "wire", denoted G = G2<Gi. Each cavity consists 
of a pair of mirrors, one of which is perfectly reflecting (shown solid) while the 
other is partially transmitting (shown unfilled). The partially transmitting 
mirror enables the light mode inside the cavity to interact with an external 
light field, such as a laser beam. The external field is separated into input 
and output components by a Faraday isolator. The optical interconnect is 
formed when light from the output of one cavity is directed into the input of 
the other, here using an additional mirror. 



Series (also called cascade) connections of quantum optical components 
were first considered in the papers [B], [3J, and certain linear feedback net- 
works were considered in [26] . Our results extend the series connection results 
in these works by including more general interfaces, and by introducing an 
efficient algebraic tool, the series product. We also introduce another prod- 
uct, which we call the concatenation product G = Gi EB G2, that is useful for 
assembling and representing systems without necessarily having connections. 
Both products may be used to describe a wide range of open quantum phys- 
ical systems (including those with physical variables that evolve nonlinearly) 
and networks of such systems (with boson field interconnects such as optical 
beams or phonon vibrations in materials). We believe our modeling frame- 
work is of fundamental system-theoretic interest. The need for general and 
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efficient methods for describing networks of quantum components has been 
recognized to some extent and has begun to emerge in the quantum optics and 
quantum information and computing literature, e.g. [27], [6], [3], [H Chap- 
ter 12], [2QI Chapter 4], [26], [I]. It is expected that an effective quantum 
network theory will assist the design of quantum technologies, just as electri- 
cal network theory and block diagram manipulations help engineers design 
filters, control systems, and many other classical electrical systems. 

Series connections provide the foundation for some important develop- 
ments in quantum feedback control, e.g. [24], [25], [23], [26], [TTJ, [12], [T3] . 
To illustrate the power and utility of our quantum network modeling method- 
ology, we analyze several examples from this literature. The series and con- 
catenation products allow us to express these quantum feedback control and 
quantum filtering examples in a simple, transparent way (there are some 
subtle technical issues in some of the examples for which we provide expla- 
nation and references). We hope this will help open up some of the quantum 
feedback control literature to control engineers, which at present is largely 
unknown outside the physics community. A number of articles and books 
are available to help readers with the background material on which the 
present paper is based. The papers [26] and [22] provide excellent introduc- 
tions to aspects of the quantum models we use. The paper [2] is a tutorial 
article written to assist control theorists and engineers by providing intro- 
ductory discussions of quantum mechanics, open quantum stochastic models, 
and quantum filtering. The book [5] is an invaluable resource for quantum 
noise models and quantum optics, while the book [21] provides a detailed 
mathematical treatment of the Hudson-Parthasarathy theory of the quan- 
tum stochastic calculus. The book [19] is a standard textbook on quantum 
mechanics. 

We begin in section [2] by discussing an analog of our results in the context 
of classical linear systems theory, elaborating further on the discussion at 
the beginning of this section. In section [3] we provide a review of some 
example quantum components (including the cavity mentioned above) and 
connections. This section includes a brief discussion of quantum mechanics, 
introduces examples of parametric representations, and provides a glimpse 
of how the general theory can be used. Open quantum stochastic models 
are described in more detail in section HI The main definitions and results 
concerning the concatenation and series products are given in section [5j in 
particular, the principle of series connections, Theorem 15.51 In general the 
series product is not commutative, but we are able to show how the order 
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can be interchanged by modifying one of the components, Theorem 15.61 A 
selection of examples from the quantum control literature are analyzed in 
section [61 The appendices contain proofs of some of the results and some 
additional technical material. 

Notation. In this paper we use matrices M = {m^} with entries rriij 
that are operators on an underlying Hilbert space. The asterisk * is used to 
indicate the Hilbert space adjoint A* of an operator A, as well as the complex 
conjugate z* = x — iy of a complex number z = x + iy (here, i = and 
x, y are real). Real and imaginary parts are denoted Re (2) = [z + z*)/2 and 
lm(z) = —%{z — z*)/2 respectively. The conjugate transpose of a matrix 
M is defined by = {m*^}. Also defined are the conjugate M" = {m*j} 
and transpose M T = {rriji} matrices, so that = (M 1 ")" = (M") T . In the 
physics literature, it is common to use the dagger f to indicate the Hilbert 
space adjoint. The commutator of two operators A, B is defined by [A, B] = 
AB — BA. S(-) is the Dirac delta function, and 5jk is the Kronecker delta. 
The tensor product of operators A, B defined on Hilbert spaces H, G is an 
operator A ® B defined on the Hilbert space H ® G (tensor product of Hilbert 
spaces) defined by (A ® B)(ip ® (j>) = (Aip) <g> (B<f)) for ip e H, e G; we 
usually follow the standard shorthand and write simply AB = A®B for the 
tensor product, and also A = A <g) I and B = I <g) B. 



2 Classical Linear Systems 

As mentioned in the Introduction (section [T]), it is common practice in clas- 
sical linear control theory to perform manipulations of block diagrams. Such 
manipulations, of course, greatly assist the analysis and design of control sys- 
tems. To assist readers in interpreting the main quantum results concerning 
series and concatenation products (section EJ), we describe concatenation and 
series products for familiar classical linear systems in algebraic terms. 
Consider two classical deterministic linear state space models 

cc j — jOu j I ^3 jU j 

y. = CjXj + DjUj (1) 

where j = 1,2. As usual, Xj, Uj and yj are vectors and Aj, Bj, Cj and Dj 
are appropriately sized matrices. These systems are often represented by the 
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matrix 



d Di 



(2) 



or the transfer function Qj(s) = Cj(sl — Aj)~ 1 Bj 

In modeling networks of such systems, one may form the concatenation 
product 
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see Figure |3j In terms of transfer functions, the concatenation of two sys- 
tems is Q(s) = diag{^!(s), Q 2 ( s )}- The concatenation product simply assem- 
bles the two components together, without making any connections between 
them. It is not a parallel connection. 




Figure 3: Concatenation product. 



Of considerable importance is the series connection, described by series 
product 

g = g 2 <Gi = ( ( b 2 c x a 2 ) ( b 2 d x ) \ , 

\ ( D 2 C, C 2 ) D 2 Di J 

see Figure [H Here the connection is specified by u 2 = yi, and so we require 
dimw 2 =dimyi. In the frequency domain, the series product is given by the 



6 



matrix transfer function product Q(s) = Q 2 ( s )Gi(s). This product describes 
a series (or cascade) connection which is fundamental to feedforward and 
feedback control. 

Notice that both products are defined in terms of system parameters 
(state space parameters or transfer function matrices). 



3 Example Components and Connections 

3.1 Some Introductory Quantum Mechanics 

Central to quantum mechanics are the notions of observables X, which are 
mathematical representations of physical quantities that can (in principle) 
be measured, and state vectors if), which summarize the status of physical 
systems and permit the calculation of expectations of observables. State 
vectors may be described mathematically as elements of a Hilbert space H, 
while observables are self-adjoint operators on H. The expected value of an 
observable X when in state if) is given by the inner product (if),Xif>). 

A basic example is that of a particle moving in a potential well, [T9| 
Chapter 14] . The position and momentum of the particle are represented by 
observables Q and P, respectively, defined by 

= #(<z), (P^){q) = -ih^{ q ) 

for if) G H = L 2 (R). Here, % = \J — 1, h = h/2n, h is Planck's constant, and 
q G R represents position values. In following subsections we use units such 
that h — 1, but retain it in our expressions in this subsection. The position 
and momentum operators satisfy the commutation relation [Q, P] = ih. The 
dynamics of the particle is given by Schro dinger's equation 

^fv(t) = HV(t), 

with initial condition V(0) = I, where H — ^~ + \muj 2 Q 2 is the Hamiltonian 
(here, m is the mass of the particle, and u is the frequency of oscillation) . The 
operator V(t) is unitary (V*(t)V(t) = V(t)V*(t) = I, where / is the identity 
operator, and the asterisk denotes Hilbert space adjoint) — it is analogous to 
the transition matrix in classical linear systems theory. State vectors and 
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observables evolve according to 



^ = V(t)i> e H, X(t) = V*(t)XV(t). 

These expressions provide two equivalent descriptions (dual), the former is 
referred to as the Schrddinger picture, while the latter is the Heisenberg 
picture. In this paper we use the Heisenberg picture, which is more closely 
related to models used in classical control theory and classical probability 
theory. In the Heisenberg picture, observables (and more generally other 
operators on H) evolve according to 

jX{t) = - l -[X{t),H(t)l (3) 

where H(t) = V*(t)HV(t). 

Energy eigenvectors ip n are defined by the equation Hip n = E n ip n for real 
numbers E n . The system has a discrete energy spectrum E n = (n + \)fajj, 
n = 0,l, 2,.... The state tp Q corresponding to E is called the ground state. 
The annihilation operator 



moj . ^ P 

{Q + i- t 



2h Imuj 

and the creation operator a* lower and raise energy levels, respectively: 
aip n = y/nip n -\, and a*ip n = \/n + l^n+i- They satisfy the canonical com- 
mutation relation [a, a*} = 1. In terms of these operators, the Hamiltonian 
can be expressed as H = ftuj(a*a + ~). Using the annihilation operator 
evolves according to 

— a(t) = —iua(t) (4) 

(JjL 

with solution a(t) = e~ luJt a. Note that also a*(t) = e luJt a*, and so commuta- 
tion relations are preserved by the unitary dynamics: [a(t), a*(t)] = [a, a*} = 
1. Because of the oscillatory nature of the dynamics, this system is often 
refereed to as the quantum harmonic oscillator. 

It can be seen that the Hamiltonian if is a key "parameter" of the quan- 
tum physical system, specifying its energy. 

3.2 Optical Cavities 

A diagram of an optical cavity is shown in Figures HI [51 together with a sim- 
plified representation. It consists of a pair of mirrors; the left one is partially 
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transmitting (shown unfilled), while the right mirror is assumed perfectly 
reflecting (shown solid). Between the mirrors a trapped electromagnetic (op- 
tical) mode is set up, whose frequency depends on the separation between the 
mirrors. This mode is described by a harmonic oscillator with annihilation 
operator a (an operator acting on a Hilbert space H (as in subsection 13. ip . 
called the initial space). The partially transmitting mirror affords the oppor- 
tunity for this mode to interact with an external free field, represented by a 
quantum stochastic process b(i) (to be discussed shortly). When the external 
field is in the vacuum state, energy initially inside the cavity mode may leak 
out, in which case the cavity system is a damped harmonic oscillator, [8]. 



input beam 



cavity 



partially reflecting 
g transmitting mirror 

mirror 

1 F output 
beam 

Figure 4: A cavity consists of a pair of mirrors, one of which is perfectly 
reflecting (shown solid) while the other is partially transmitting (shown un- 
filled). The partially transmitting mirror enables the light mode inside the 
cavity to interact with an external light field, such as a laser beam. The 
external field is separated into input and output components by a Faraday 
isolator. 



Quantization of a (free) electromagnetic field leads to an expression for 
the vector potential 

A(x,t) = J K(w)[6(w)e-^ +ira/c + &*(w)e !wi -^ /c ](iw, 

for a suitable coefficients k(u), and annihilation operators b(u). Such a field 
can be considered as an infinite collection of harmonic oscillators, satisfying 
the singular canonical commutation relations 

[b(u),b*(u')} = 5(u-u'), 

where 5 is the Dirac delta function. 
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Figure 5: A simplified representation of the cavity from Figure H] which omits 
the Faraday isolator. It shows input B and output B fields and the cavity 
mode annihilation operator a. This representation will be used for the re- 
mainder of this paper. 



An optical signal, such as a laser beam, is a free field with frequency con- 
tent concentrated at a very high frequency ujq ~ 10 14 rad/sec. The fluctua- 
tions about this nominal frequency can be considered as a quantum stochastic 
process consisting of signal plus noise, where the noise is of high bandwidth 
relative to the signal. Indeed, a coherent field is a good, approximate, model 
of a laser beam, and can be considered as the sum b(t) = s(t) + bo{t), where 
s(t) is a signal, and b (t) is quantum (vacuum) noise. Such "signal plus 
noise" models are of course common in engineering. 

The cavity mode-free field system has a natural input-output structure, 
where the free field is decomposed as a superposition of right and left travel- 
ing fields. The right traveling field component is regarded as the input, while 
the left traveling component is an output, containing information about the 
cavity mode after interaction. The interaction facilitated by the partially 
transmitting mirror provides a boundary condition for the fields. The two 
components can be separated in the laboratory using a Faraday isolator. This 
leads to idealized models based on rotating wave and Markovian approxima- 
tions, where, in the time domain, the input optical field (when in the ground 
or vacuum state) is described by quantum white noise b(t) = b (t) 0, Chap- 
ters 5 and 11], which satisfies the singular canonical commutation relations 



In order to accommodate such singular processes, rigorous white noise and 
Ito frameworks have been developed, where in the Ito theory one uses the 




(5) 
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integrated noise, informally written 



B(t) = [ b(s)ds 



The operators B(t) are defined on a particular Hilbert space called a Fock 
space, F, j21j sec - 19]- When the field is in the vacuum (or ground) state, 
this is the quantum Wiener process which satisfies the Ito rule 



(all other Ito products are zero). Field quadratures, such as B(t) + B*(t) and 
—i(B(t) — B*(t)) are each equivalent to classical Wiener processes, but do not 
commute. A field quadrature can be measured using homodyne detection, [8j 
Chapter 8]. 

The cavity mode-free field system can be described by the Hamiltonian 



where the first term represents the self-energy of the cavity mode (the number 
A is called the "detuning" , and represents the difference between the nominal 
external field frequency and the cavity mode frequency), while the remaining 
two terms describe the energy flow between the cavity mode and the free 
field (a photon in the free field may be created by a loss of a photon from the 
cavity mode, and vice versa). This Hamiltonian is defined on the composite 
Hilbert space, the tensor product H ® F; the tensor product is not written 
explicitly in the expression (jH]). 

The Schrodinger equation for the cavity-free field system is derived from 
dBl) under certain assumptions [8], and is given by the Ito quantum stochastic 
differential equation (QSDE) 



with vacuum input and initial condition V(0) = I, so that V(t) is unitary. 
The complete cavity mode-free field system thus has a unitary model. In the 
Heisenberg picture, cavity mode operators X (operators on the initial space 



dB(t)dB*(t) = dt 





dV(t) = {y/^adB*(t) - y/ja*dB(t) 
——a*adt — iAa*adt}V(t), 



(7) 
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H) evolve according the quantum Ito equation 



dX(t) = -iA[X(t),a*(t)a(t)]dt (8) 
+ l(a*{t)[X(t),a(t)] + [a*{t),X(t)]a(t))dt 
+ yfydB*(t)[X(t),a(t)] + yfy[a*(t),X(t)]dB(t). 

Here, 7 > is a parameter specifying the coupling strength, and is related 
an approximation of the function k(u) in the Hamiltonian (jB]). In particular, 
for X = a, the cavity mode annihilation operator, we have 

da{t) = -(^ + iA)a(t)dt-^(dB(t); (9) 

cf. (|1J). The output field B(t) is given by 

dB(t) = y/ja(t)dt + dB(t), (10) 

where one can see the "signal plus noise" form of the field. 

This is an example of an open quantum system, characterized by the 
parameters and Aa*a; the latter being the cavity mode Hamiltonian 
(specifying internal energy), and the former being the operator coupling the 
cavity mode to the external field (specifying the interface). These parameters 
are operators defined on the initial space H. These parameters specify a 
simpler, idealized model employing quantum noise, in place of the more basic 
but complicated Hamiltonian (jHJ). 



3.3 Optical Beamsplitters 

A beamsplitter is a device that effects the interference of incoming optical 
fields Ai,A 2 and produces outgoing optical fields Ai,A 2 , Figure El The 
relationship between these fields is 

A 1 {t)=/3A 1 (t)-aA 2 (t), A 2 {t)=aA 1 (t)+pA 2 (t), (11) 

where a and (3 are complex numbers describing the beamsplitter relations, 
and they satisfy a* a + (3* (3 = 1, a* (3 = a/3* (here the asterisk indicates the 
conjugate of a complex number). 

The initial space is trivial, H = C, the complex numbers; nevertheless, 
the Schrodinger equation for the beamsplitter is 

dV{t) = {(S-I)dA}V{t), (12) 
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Figure 6: Diagram of an optical beamsplitter showing inputs Ai,A% and 
outputs Ai,A 2 fields. 



with initial condition V(0) = /, where S is the unitary matrix defined by 
( THj) below, I is the identity matrix, and A is the matrix of gauge processes 

Here, A^ describes the destruction of a photon in channel j and the cre- 
ation of a photon in channel i. In terms of their formal derivatives, Aij(t) = 
a*(s)dj(s)ds, where Ai{t) = J^ai(s)ds. The self-adjoint processes Ajj are 
equivalent to classical Poisson processes when the channels are in coherent 
states (signal plus quantum noise). These counting processes may be ob- 
served by a photodetector, [8], Chapters 8 and 1 1] . 

This open system is characterized by the unitary parameter matrix 

s =(f~;)- (i4) 

which describes scattering among the field channels. The matrix S specifies 
the interface for the beamsplitter. 

3.4 Open Quantum Systems 

In general, as we shall explain in more detail in section [4j open quantum 
systems with multiple field channels are characterized by the parameter list 

G=(S,L,if) (15) 

where S is a square matrix with operator entries such that S^S = SS^ = I 
(recall the notational conventions mentioned at the end of section [T]), L is a 
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column vector with operator entries, and if is a self-adjoint operator. The 
matrix S is called a scattering matrix, the vector L is a coupling vector, 
together, these parameters specify the interface between the system and the 
fields. The parameter H is the Hamiltonian describing the self-energy of the 
system. Thus the parameters describe the system by specifying energies — 
internal energy, and energy exchanged with the fields. All operators in the 
parameter list are defined on the initial Hilbert space H for the system. 

The closed, undamped, harmonic oscillator of subsection 13.11 is specified 
by the parameters 

H = (_,_,cja*a) (16) 

(the blanks _ indicate the absence of field channels), while the open, damped 
oscillator (cavity) of subsection 13.21 has parameters 

C = (I, y/^a, Aa*a). (17) 

The beamsplitter, described in subsection 13.31 has parameters 

M = ((f -;),0,0). (18) 
3.5 Series Connection Example 

B 2 = B 2 = A 2 




a 



Figure 7: Beam splitter (left) and cavity (right) network. 

Consider the feedforward network shown in Figure [7J where one of the 
beamsplitter output beams is fed into an optical cavity. From the previ- 
ous subsections, we see that the quantum stochastic differential equations 
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describing the network are 



da(t) 


= (-- + iA)a(t)dt - yfi dB x {t) 


Mt) 


= pAx{t) - aA 2 (t) 


Mt) 


= aA 1 (t)+ (3A 2 {t) 


Bi(t) 


= Mt) 


B 2 (t) 


= A 2 (t) 


dB^t) 


= ^/ja(t)dt + dBt(t) 


dB 2 (t) 


= dB 2 (t). 



(19) 

(20) 
(21) 
(22) 
(23) 
(24) 
(25) 

It can be seen that algebraic manipulations are required to describe the 
complete system (in general such manipulations may be simple in principle, 
but complicated in practice). The key motivation for this paper is more 
efficient algebraic methods for describing such networks. 

We now describe how the parameters for the complete network may be 
obtained. We first assemble the field channels into vectors as follows: 



A 2 ) ' V B 2 J ' V A 2 J ' \B 2 

The beamsplitter acts on the input vector A, and is described by the pa- 
rameters M given in equation (fT8{) ). Now the beamsplitter output has two 
channels, while the cavity has one channel (described by the parameters C, 
equation (TT7I) ). and so we augment the cavity to accept a second channel in 
a trivial way. This is achieved by forming the concatenation C ffl N, where 
N = (1,0,0) represents a trivial component (pass-through). The augmented 
cavity C ffl N can now accept the output of the beamsplitter, so that the 
complete network is described as a series connection as follows: 

G=(CfflN)«M. (26) 

The definition of the concatenation ffl and series < products will be explained 
below in section 151 (Definitions 15.11 and 15.31 and the principle of series connec- 
tions, Theorem I5.5p . By applying these definitions, we obtain the network 
parameters 

G=(('l - a ),(^ a ),Aa'a). (27) 



a (3 J ' V 
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B 2 



Figure 8: Beam splitter-cavity network representation illustrating the net- 
work model given by (1271) . 



A schematic representation of the network is shown in Figure [H which illus- 
trates the important point that components, parts of components, as well as 
the complete network, are described by parameters of the form [T5\) . 

For the purposes of network modeling and design, it can be useful to 
perform manipulations of the network to yield equivalent networks; this, of 
course, is common practice in classical electrical circuit theory and control 
engineering. For instance, in our example we could move the beam split- 
ter to the output, but the cavity should be modified (to have two partially 
transmitting mirrors) as follows (see Remark 15. 7\) : 

G = (CfflN)<M = M<(C / fflN / ). (28) 

Here, the modified cavity C EB N' (see Figure [H]) is described by the subsys- 
tems 

C' = {I,p*y/ja, Aa*a) , N' = (I, -a'/j-a, 0) . (29) 



-Bl = A X 


modified cavity 


Bl 

An ' 


y b 2 















Figure 9: Equivalent beam splitter and cavity network. 



The connections described here so far are unidirectional field mediated 
connections. Components interact indirectly via a quantum field, which acts 
as a quantum "wire" . One can also consider bidirectional direct connections, 
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which can be accommodated by using interaction Hamiltonian terms in the 
models. Our emphasis in this paper will be on field mediated connections, 
with direct connections readily available in the modeling framework if re- 
quired. See subsection 15.41 



4 Open Quantum Stochastic Models 



In this section we describe in more detail the open quantum models of the 
type encountered in section [3j Specifically, we consider models specified by 
the parameters G = (S,L, H) (recall (fl5]l ). where 



Sin ' 



s, 



( L, \ 



are respectively a scattering matrix with operator entries satisfying S^S = 
SS^ = I, and coupling vector with operator entries, and if is a self-adjoint 
operator called the Hamiltonian (this parameterization is due to Hudson- 
Parthasarathy, [15] , and is closely related to a standard form of the Lindblad 
generator, given in (133]) below). The operators constituting these parameters 
are assumed to be defined on an underlying Hilbert space H, called the initial 
space. These parameters specify an open quantum system coupled to n field 
channels with corresponding gauge processes: 
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All differentials shall be understood in the Ito sense - that is, dX (t) = 
X (t + dt) — X (t). We assume that these processes are canonical, mean- 
ing that we have the following non-vanishing second order Ito products: 
dAj (t) dA k (t)* = S jk dt, dA jk (t)dAi(t)* = S u dAj(t)*, dAj (t) dA kl (t) = 
S jk dAi(t) and dA jk (t) dA tm (t) = 6 kl dA jm (t). 

If we consider the open system specified by G = (S, L, H) with canonical 
inputs, the Schrodinger equation 



dV(t) = {tr[(S - I)dA] + dA f L 



(30) 



- VSdA 



-VLdt 

2 



iHdt}V{t) = dG(t) V(t) 
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with initial condition V(0) = I determines the unitary motion of the system. 
Equation (1311 serves as the definition of the time-dependent generator dG(t). 
Given an operator X defined on the initial space H, its Heisenberg evolution 
is defined by 

X(t)=i t (X)=V(t)*XV(t) (31) 

and satisfies 

dX(t) = (C m (X(t)) - i[X(t),H(t)])dt 
+c/A t (t)S t (t)[X(t),L(t)] + [V(t),X(t)]S(t)dA(t) 

+tr[(S f (*)X(t)S(t) - X(t))dA(t)]. (32) 

In this expression, all operators evolve unitarily according to (l3"Tj) (e.g. L(t) = 
j 4 (L)) (commutators of vectors and matrices of operators are defined component- 
wise), and tr denotes the trace of a matrix. We also employ the notation 

£ L (X) = ±V[X,L] + ±[V,X]L 

n 

= E(2^ [X '^' ] + 2 [ ^' X]Lj); (33) 
i=i 

this is called the Lindblad superoperator in the physics literature (it is analo- 
gous to the transition matrix for a classical Markov chain, or the generator of 
a classical diffusion process). The dynamics is unitary, and hence preserves 
commutation relations. The output fields are defined by 

A(i) = V*{t)A(t)V{t), A(t) = V*(t)A{t)V(t), (34) 

and satisfy the quantum stochastic differential equations 

dk(t) = S(t)dA(t)+L(t)dt 
dA(t) = S tt (t)dA(t)S T (t) + S t {t)dA t (t)L T {t) 
+L i {t)dA T {t)S T {t) + l)(t)l/(t)dt, 

where L(t) = jt(L), etc, as above. The output processes also have canonical 
quantum ltd products. 

In the physics literature, it is common practice to describe open sys- 
tems using a master equation (analogous to the Kolmogorov equation for the 
density of a classical diffusion process) for a density operator p, a convex 
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combination of outer products ipip* (here ip is a state vector). Master equa- 
tions can easily be obtained from the parameters G = (S, L, H); indeed, we 
have 

j t P = i[p,H(t)]+£' m (p), (35) 

where C' L (p) = L T pL tl — |L"L T p — ipL"L T is the adjoint of the Lindbladian: 
ti[p{t) C\,{X)] = tr[£' L (p) X}. Note that while the master equation does not 
depend on the scattering matrix S, this matrix plays an important role in 
describing the architecture of the input channels, as in subsections 13.51 and 
I6.2L We also mention that if an observable of one or more output channels 
is continuously monitored, then a quantum filter (also called a stochastic 
master equation) for the conditional density operator can be written down 
in terms of the parameters G = (S, L, H); an example of this is discussed in 
subsection 16. 3[ see [2]. 

Open systems specified by parameters G = (S, L, H) preserve the canon- 
ical nature of the quantum signals. However, if the inputs are not canonical, 
one will need to modify the equations for the unitary, the Heisenberg dy- 
namics, and the outputs, etc, to accommodate non-canonical correlations; 
we do not pursue this matter further here, and in this paper we will always 
use canonical quantum signals. 

5 The Concatenation and Series Products and 
their Application to Quantum Networks 

This section contains the main results of the paper. The concatenation and 
series products are defined in subsection 15. 1[ and applied to a feedback ar- 
rangement in Theorem 15.51 the principle of series connections (subsection 
15.21) . This is followed in subsection 15.31 with a specialization to cascade net- 
works, and a consideration in subsection 15.41 of reducible networks. These 
results are applied to a range of examples in section 

5.1 Definitions 

In this subsection we define two products between system parameters. It is 
assumed that both systems are defined on the same underlying initial Hilbert 
space, enlarging if necessary by using a tensor product. 
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Definition 5.1 ( Concatenation product) Given two systems Gi = (Si, Li, Hi) 
and G 2 = (S 2 ,L 2 ,i? 2 ), we define their concatenation to be the system Gi EB 
G 2 by 

GifflG 2 = (( S - + (36) 

The concatenation product is useful for combining distinct systems, or 
for decomposing a given system into subsystems. It does not describe in- 
terconnections via field channels, but does allow for direct connections via 
the Hamiltonian parameters. Systems without field channels are included by 
employing blanks; set (_, _, H) EB (_, _, H') := (_, _, H + H') and more generally 
(_, _, H) BE! (£", I/, H') = (S', L', H>) EB (_, _, H) := (S", L', # + #'). 

Definition 5.2 (Reducible system) We say that a system G = (S, L, H) is 
reducible if it can be expressed as 

G = Gi EB G 2 (37) 

for two systems G\ and G 2 . In particular, the parameters of a reducible 
system have the form 

S =(o S 2 )■ L =(l0' n=H l + H,. (38) 

Such decompositions are not unique. Furthermore, if one or more of the 
subsystems is reducible, the reduction process may be iterated to obtain a 
decomposition G = EB^Gj. 

Definition 5.3 (Series product) Given two systems G\ = (Si,Lq,iJi) and 
G 2 = (S 2 , L 2 , H 2 ) with the same number of field channels, the series product 
G 2 < Gi defined by 

G 2 < Gi = (S 2 Si, L 2 + S 2 Li, 

H 1 + H 2 + ^(I4S 2 L! - L{S 2 L 2 )). 

As will be explained in the following subsection, the series product speci- 
fies the parameters for a system formed by feeding the output channel of the 
first system into the input channel of the second. Both of these products are 
powerful tools for describing quantum networks. 
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Remark 5.4 Let dGj(t) denote the infinitesimal Ito generators correspond- 
ing to parameters Gj = (S^, Ly, Hj), for j — 1,2 respectively, as constructed 
in ()30p . The generator corresponding to G2 < Gi is then 

dG{t) = dGx{t) + dG 2 {t) + dG 2 (t)dG l {t). (39) 

The last term is to be computed using the Ito table for second order products 
of differentials. 



5.2 Feedback 

Let us consider a reducible system G = Gi EH G 2 (recall Definition 15. 2p . 
where number of channels in the factors is the same (i.e. dim hi = dim L 2 ). 
The setup is sketched in Figure [10J We investigate what will happen if we 
feed one of the outputs, say A x back in as the input A 2 . Either of the two 
diagrams in Figure [TT] may serve to describe the resulting feedback system. 
Note that the outputs will be different after the feedback connection has been 
made. 

















— 2~ 





Figure 10: Reducible system GiEHG 2 with inputs A 1; A 2 and outputs A 1; A 2 . 



1 — ► 
























A 2 = Ai 



A 2 



At. 



A 2 



Figure 11: Direct feedback system G 2 < Gi, with input Ai and output A 2 . 



We now state our main result applying the series product to feedback. 
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Theorem 5.5 (Principle of Series Connections) The parameters G 2 ^i 
for the feedback system obtained from G± EH G 2 when the output of the first 
subsystem is fed into the input of the second is given by the series product 
G 2 ^i = G 2 < Gi. 

A proof of this theorem is given in Appendix [Bl 
5.3 Cascade 

In our treatment of series connections, we nowhere assumed that the ma- 
trix entries commuted, and this of course facilitated feedback. However, the 
principle of series connections also applies to the special case where the sub- 
systems commute, as in a cascade of independent systems, as shown in Figure 

E2J0 

To formulate the cascade arrangement, we first consider the concatenation 
of the two systems Gi EH G 2 . The system G = Gi EH G 2 is reducible with 
components Gj. 



A, 




Ai = A 


2 


A 




— • — 




— • — 





Gi G 2 



Figure 12: Cascade of independent quantum components, G 2 < Gi. 

The notion of cascaded quantum systems goes back to Carmichael [3], 
who used a quantum trajectory analysis, and Gardiner [6] who used (scalar) 
quantum noise models of the form Gj = (l,Lj,Hj) (no scattering). As 
a special case of the series principle, we see that the cascaded generator for 
this type of setup is G cascade = G 2 «Gi = (1, L 1 + L 2 , H 1 + H 2 + Im {Lj^Li}). 
This is entirely in agreement with Gardiner's analysis, cf. [HJ Chapter 12] with 
Lj = y^yjcj where we have L 2 <~± = y/jiCi + v /t 2 c and H 2 ^i = H\ + H 2 + 

^v/7i72"(c£ci -c*c 2 ). 

We now consider cascade arrangements and ask what happens if we try to 
swap the order of the components. Since the series product is not in general 

1 Indeed, the reason we use the term "series" is to indicate that it applies more generally 
than to cascades of independent components. 
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commutative, we cannot expect to be able to swap the order without, say, 
modifying one of the components. We now make this precise as follows. 

We say that two systems are parametrically equivalent if their parameters 
are identical. This implies that, for the same input, they produce the same 
internal dynamics and output. Consider the cascaded systems shown in 
Figure E3 



G 5 



G 



G 1 



Figure 13: Equivalent Systems. 



We assume that the initial inputs are canonical in both cases and ask, for 
fixed choices of Gi and G 2 , what we should take for G 2 so that the setups 
are parametrically equivalent. 

Theorem 5.6 The two cascaded systems shown in Figure [731 are parametri- 
cally equivalent if and only if 

G 2 <G 1 = G 1 < G' 2 . (40) 

Furthermore, if (Sj, Lj, Hj) are the parameters for Gi and G 2 (j = 1,2) , 
then the parameters (S' 2 ,L 2 , H' 2 ) of G' 2 are uniquely determined by 

S' 2 = SjS2Si, 

L 2 = S| (S2 — I) Li + S|Tj2, 

H' 2 = iJ 2 + Im|L5(S2 + I)L 1 -L t 1 S 2 L 1 }. (41) 
The proof of this theorem is given in Appendix 

Remark 5.7 A useful special case of this result is moving a scattering matrix 
from the input to the output of a modified system: 

(S, L, H) = (I, L, H) < (S, 0, 0) = (S, 0, 0) < (I, S f L, H). (42) 

This is illustrated in subsection 13.51 □ 
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5.4 Reducible Networks 



Networks can be formed by combining components with the concatenation 
and series products. Within this framework, components may interact di- 
rectly, or indirectly via fields. This framework is useful for modeling existing 
systems, as we have seen above, as well as for designing new systems. 

Let {Gj} be a collection of components, which we may combine together 
to form an unconnected system G = ESjGj. The components may interact 
directly via bidirectional exchanges of energy, and this may be specified by 
a direct connection Hamiltonian K of the form 

K = tJ2(N* k M k -M* k N k ), (43) 

k 

where M k , N k are operators defined on the initial Hilbert space for G. The 
components may also interact via field interconnects, specified by a list of 
series connections 

y = {G h «G kl ,...,G jn «G kn } (44) 

such that (i) the field dimensions of the members of each pair are the same, 
and (ii) each input and each output (relative to the decomposition G = 
BBjGj) has at most one connection. 

A reducible network N is the system formed from G by implementing 
the connections (|43|) and (JHJ). The parameters of the network N may be 
obtained as follows. A series chain is a system of the form 

C = G jl <G h <---<G jm <G km . 

Let ^ denote the set of maximal-length chains drawn from the list of series 
connections (JHJ), and let % denote the set of components not involved in 
any series connection. Then the reducible network is given by 

N = (ffl Gfc6 *G fc ) ffl {m Ci evCj) ffl (1, 0, K). (45) 

An example of a reducible network is shown in Figure [TH 

Remark 5.8 The examples considered in section below are all important 
examples of reducible networks that have appeared in the literature. How- 
ever, we mention that there are important examples of quantum networks 
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G N 

Figure 14: A reducible network N = Gi EB (G4 < G3 < G2) formed from the 
collection G = Gi EB G2 EB G3 EB G4 of components with connections specified 
by the list of series connections 5? = {G 3 < G 2 , G 4 < G 3 }. 



that are not reducible. An example of a non-reducible network was consid- 
ered by Yanagisawa and Kimura, [261 Fig- 4], which consists of two systems 
in a feedback arrangement formed by a beam splitter, as occurs if in Figure 
U\ we connect the output B\ to the input A 2 (i.e. setting A 2 = B\). The 
feedback loop formed in this way is "algebraic", and the resulting in-loop 
field is not a free field in general. A general theory of quantum feedback 
networks, both reducible and non-reducible, is given in [11]. □ 

6 Examples 

In this section we look at a number of examples from the literature which 
can be represented by reducible networks. 

6.1 All-Optical Feedback 

We consider a simple situation first introduced by Wiseman and Milburn as 
an example of all-optical feedback, [251 section II. B. A]. Referring to Figure 
[T5l vacuum light field A\ is reflected off mirror 1 to yield an output beam A\ 
which results from interaction with the internal cavity mode a. This beam 
is reflected onto mirror 2, as shown, where it constitutes the input A 2 . It is 
assume that both mirrors have the same transmittivity, so that we can model 
the coupling operators for the two field channels as L\ = L 2 = w^ya, where 

7 is the damping rate. We may also assume that the light picks up a phase 
S = e ld when reflected by the cavity mirror. 
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cavity 



A 2 



mirror 1 



mirror 2 



optical interconnect 



Figure 15: All-optical feedback for a cavity. The feedback path is a light 
beam from mirror 1 to mirror 2, both of which are partially transmitting). 
There is a phase shift 9 along the feedback path. 



Before feedback, the cavity is described by 



G = (I, I ^ J,0) = (l,L 1 ,0)ffl(l,L 2 ,0). 
The phase shift between the mirrors is described by the system (S, 0, 0). 



Mi 



An 



1 






A 2 










-~v 2 




(5,0,0) 



(5,0,0) 

Figure 16: Representations of the all-optical feedback scheme of Figure [TS] 
as reducible networks. 



Two equivalent reducible network representations are shown in Figure [ 
From the left diagram in Figure [TU we see that the closed loop system is 
described by 

G d = (1,L 2 ,0)<(5,0,0)<(1,L 1 ,0) 

= (5,5L 1 + L 2 ,^(L^L 1 -L^*L 2 )). 
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Here we have twice applied the formulas (1391) given in Definition 15.31 

Alternatively, we may use our theory of equivalent components (Theorem 
15.61) to move the phase change (S, 0, 0) to the very end, as shown in the right 
diagram in Figure [161 Then 

G cl = (S,0,0)<(1,S*L 2 ,0)<(1,L 1 ,0) 

= {S,SLx + L 2 ^{L* 2 SLx-L{S*L 2 )\ 

as before. Either way, the closed loop feedback system is described by G c i = 
(S ch L ch H ci ) where 

S d = S = e ie , 

L cl = SLx + L 2 = (l+e 4e ) v^a, 
H c \ = Im {L^SLi} = 7 sin 9 a^a. 

From this we obtain the Heisenberg dynamical equation for the cavity mode 

da = - [a, (1 + e ie ) ^/ja^] dA l 

~ (l + e ld ) (l + e~ ie ) adt - z 7 sin 9 adt 
= - (1 + e ie ) {^jdA 1 + jadt) , 

and the input/output relation, in agreement with (25J eq. (2.29)], 

dA 2 = e ie dA 1 + (1 + e ie ) yfyadt. 

6.2 Direct Measurement Feedback 

In the paper |24j . Wiseman considers two types of measurement feedback, one 
involving photon counting, and another based on quadrature measurement 
using homodyne detection (which is a diffusive limit of photon counts). In 
both cases proportional feedback involving an electrical current was used. 
We describe these feedback situations in the following subsections using our 
network theory. 

Consider the measurement feedback arrangement shown in Figure [XT], 
which shows a vacuum input field A, a control signal c, a photodetector PD, 
and a proportional feedback gain k. 
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Figure 17: Direct feedback of photocurrent obtained by photon counting 
using a photodetector (PD). 



Before feedback, the quantum system is described by 

G = (1,L,H + Fc), (46) 

where Hq and F are self-adjoint, and c represent a classical control variable. 
The photocurrent j(t) resulting from ideal photodetection of the output field 
is given by 

u j(t)dt" =dA + LdA^ + tfdA + tfLdt, (47) 

where, mathematically, the photocurrent j(t) is the formal derivative of a 
field observable (a self-adjoint commutative jump stochastic process) A(t) 
(the output gauge process) whose Ito differential is given by the RHS of 
( )47l) . The feedback is given by 

c(t)=kj(t), (48) 

where k is a (real, scalar) proportional gain. The feedback gain can be 
absorbed into F, and so we assume k = 1 in what follows. 

An alternative is to again consider the quantum system G given by (T46l) , 
but replace the photodetector PD in Figure [T7] with a homodyne detector 
HD@ The homodyne detector then produces a photocurrent j(t) given by 

u j(t)dt" = dJ(t) = (L(t) + l)(t))dt + dA(t) + dA\t). 

2 An ideal homodyne detector HD takes an input field A and produces a quadrature, 
say A + A* (real quadrature), thus effecting a measurement. This is achieved routinely to 
good accuracy in optics laboratories, [51 Chapter 8]. 
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The feedback is given by (j48p as above, with feedback gain absorbed into 
F, as above. The measurement result J(t) is a field observable (here a self- 
adjoint commutative diffusive process). 

In order to describe these types of direct measurement feedback within 
our framework, we view the setup before feedback as being described by 



Here, Go describes the internal energy of the system and its coupling to 
the input field A. The second term, Gf b , describes the way in which the 
classical input signal is determined from a second quantum input field (which 
will be replaced by the output A when the feedback loop is closed). The 
idea is that by appropriate choice of the coupling operator Lf b , the relevant 
observable of the field can be selected. In this way, the photodection and 
homodyne detection measurements are accommodated. The singular nature 
of the feedback signal (which contains white noise in the homodyne case) 
means that care must be taken to describe it correctly. The correct form 
of the parameters is given by the Holevo parameterization (Appendix [Aj 
equation (|55l) ) rather than the expression arising from the implicit-explicit 
formalism of [21] , since the later does not capture correctly gauge couplings, 
see Appendix |A] We shall interpret the feedback interaction as being due to a 
Holevo generator K fb (t) = H 00 t-\-H iA(t)+H w A*(t)+H n A(t), see Appendix 
IA~t equation (l54"j) . The closed loop system after feedback is given by the series 
connection G d = G fb < G = (S fb , L fb + S fb L, H + H fb + Im (L* fb S fb L )) . 

6.2.1 Photon Counting 

Here we take Kf b (t) = FA(t), so that Sf b = e~ tF , see Appendix IA~| equa- 
tion flo^]) . Note that this coupling picks out the required photon number 
observable of the field. We then have Gf b = (e~ tF , 0, 0) and so 



This is illustrated in Figure [T8j The resulting Heisenberg equation agrees 
with the results obtained by Wiseman, [2H eq. (3.44)], which we write in 
our notation as 



G = ( 



1, L, Ho) EB (Sf b , Lf b , Hf b ) = Go EB G/&. 



G cl = (e-* F ,e- iF L,H ). 



dX 



(-i[X, H ] + C e -iF L (X))dt + (e iF Xe~ iF 
+e iF [X,e- iF L]dA* + [L*e iF , X}e~ iF dA. 



X)dA 



(49) 
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(Technical aside. Note that if we set E(t) = EA(t), with E self-adjoint, 
then the Stratonovich equation dV(t) = —idE(t) o V(t) = —idE(t)V(t) — 

\dE(t)dV(t) is equivalent to dV{t) = S fb dA(t)V(t) where S fb = ^f§. 

Therefore the implicit form [21] is not the Stratonovich form |10j.) 



L 



-B- 



Figure 18: Representation of the direct photocount feedback scheme of Figure 
[T71 as a reducible network. 



6.2.2 Quadrature Measurement 

Here we take K fb (t) = F(A*(t) + A(t)) in which case G fb = (1, — iF, 0), see 
Appendix El equation ( 1551) . The skew-symmetry of — iF ensures that the 
coupling selects the desired field quadrature observable. After feedback, the 
closed loop system is 

G d = (l,L-iF,H + ~(FL + L*F)) 

using (|39l) . This is illustrated in Figure [T9l The resulting Heisenberg equation 
then agrees with [!Ml eq. (4.21)], which we write as 

dX = (-t[X,H + ^(FL + L*F)}+C L ^ tF (X))dt 

+[X, (L - iF)}dA* + [(L - iF)*, X})dA. (50) 

(Technical aside. Note that for diffusions (that is, no gauge terms) the 
Holevo generator and Stratonovich generator coincide: that is, dV{t) = 

( e -idK fb (t) _ ig the 

same as dV(t) = -idK fb (t) o V(t), Appendix [Al) 
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Figure 19: Representation of the direct homodyne feedback scheme (Figure 
[IT] with HD replacing PD) reducible network. 



6.3 Realistic Detection 

Consider a quantum system G q continuously monitored by observing the 
real quadrature B + B* of an output field B. This measurement can ideally 
be carried out by homodyne detection, but due to finite bandwidth of the 
electronics and electrical noise, this measurement could be more accurately 
modeled by introducing a classical system (low pass filter) and additive noise 
as shown in Figure [201 as analyzed in [23]. Here, B is a vacuum field, I is the 
output of the ideal homodyne detector (HD), v is a standard Wiener process, 
and Y is the (integral of) the electric current providing the measurement 
information. 

We wish to derive a filter to estimate quantum system variables X q from 
the information available in the measurement Y. 




dc lcc lion system Gc 



Figure 20: Model of a realistic detection scheme for a quantum system, 
showing ideal homodyne detection followed by a classical system (e.g. low 
pass filter) and additive classical noise. 



The quantum system is given by G q = (1, L q , H q ), and the classical de- 
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tection system is given by the classical stochastic equations 

dx{t) = f(x(t))dt + g(x(t))dw(t), 
dY(t) = h(x(t))dt + dv{t), (51) 

where x(t) £ M n , y(t) £ R, f,g are smooth vector fields, h is a smooth 
real- valued function, and w and v are independent standard classical Wiener 
processes. As described in the Appendix [Dj this classical system is equivalent 
to a commutative subsystem of G c = (1, L c i, H c ) EB (l,L c2 ,0), where L cl = 
—ig T p — |V T g, L c2 = \h and H c = \(f T p + p T f)- We represent the system 
of Figure [20] as a redicible network, as shown in Figure [2TJ 
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Figure 21: Representation of the realistic detection scheme of Figure [20] as a 
reducible network. 



Here, the classical noises are represented as real quadratures w = Ai + A\, 
v = A2 + A2. Note that since L c \ is skew-symmetric, only the real quadrature 
w = A\ + A\ = B + B* affects the classical system (this captures the ideal 
homodyne detection). The complete cascade system is 

G = ((l,L cl ,# c )«(l,L 9 ,# g ))ffl(l,L c2 ,0) (52) 

= (I,( Ll 2 c2 Lcl ),H q + H e + ^(L* cl L q - L* q L cl )) 

Applying quantum filtering pQ, [2], the unnormalized quantum filter for 
the cascade system G is 

da t (X) = a t (-i[X, H q + H c + ^{L* cl L q - L* q L cl )} 

+£/ r T x (X))dt + a t (L* c2 X + XL c2 )dy. (53) 
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Here, X is any operator defined on the quantum-classical cascade system. 
For instance, X = X q <g> tp, where (p is a smooth real valued function on R n . 
In particular, if X = X q is a quantum system operator, one can compute the 
desired estimate of X q from irt(X q ) = a t (X q )/a t (l). 

Equation (133]) can be normalized, and compared with [23], eq. (17)]. In 
the case that the quantum system is a linear gaussian system, and the filter 
is a linear system, the complete filter reduces to a Kalman filter from which 
the desired quantum system variables can be estimated. 

7 Conclusion 

In this paper we have presented algebraic tools for modeling quantum net- 
works. The tools include a parametric representation for open quantum sys- 
tems, and the concatenation and series products. The concatenation prod- 
uct allows us to form a larger system from components, without necessarily 
including connections. The series product, through the principle of series 
connections (Theorem 15.51) . provides a mechanism for combining systems via 
field mediated connections. We demonstrated how to model a class of quan- 
tum networks, called reducible networks, using our theory and we illustrated 
our results by examining some examples from the literature. 

Future work will involve further development of the network theory de- 
scribed here, and applying the theory to develop control engineering tools 
and to applications in quantum technology, e.g. [T6j . 

A Time-Ordered Exponentials in the sense of 
Holevo 

Holevo [H] developed a parameterization of open system dynamics that is 
different to the Hudson-Parthasarathy parameters G = (S,L,H). Holevo's 
parameterization is defined as follows. Let 

K (t) = H 00 t + H 01 A{t) + H w A*(t) + H n A(t), (54) 

where {H a/3 } consists of bounded operators with H a/3 = Hp a , and the indices 
a, (3 range from to 1 (here we are considering a single field channel for 
simplicity). The time- ordered exponential with Holevo generator {H a p} is 
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the unitary adapted process U satisfying the quantum stochastic differential 
equation 

dU it) = (e~ ldKW -1)U (t) 
with U (0) = 1, [H], [9]. Expanding the differential e - jdX (*) — 1 we obtain 



/ — ' n! 



n>l 



Now for a system with parameters G = (S, L, H) we have 

dU(t) = {(S- I)dA(t) + LdA*(t) - L*dA 
-(iH + -L*L)dt}U{t). 

Comparing these expressions, we find that 

S = exp (-iH n ) , L = H 10 , 

-nil 

H = H 00 — H i — 11 , — U ^ Hi . (55) 

The relationship between the generating coefficients H a p and the param- 
eters G = (S, L, H) are exactly as occur in the implicit-explicit formalism 
of [23], however, this formalism only coincides with the Stratonovich-Ito cor- 
respondence in the case where Hu = [TO] . 



B Proof of Theorem 15.5 



There are a number of independent derivations of the series product. For 
instance it can be derived from a purely Hamiltonian formalism for quantum 
networks [11], alternatively Gardiner's arguments in the Heisenberg picture 
can be extended to include the scattering terms [12]. Here we present a 
discretization argument for the input/output fields based on [9]. Rather 
than considering a continuous noise source, we take a beam consisting of 
qubits (spin one-half particles) with a rate of one qubit every r seconds. 
A qubit has the Hilbert space H = C 2 spanned by a pair of orthogonal 
vectors eo and e\. We define raising/lowering operators a ± for each qubit by 
cr + (aeo+j3ei) = at\ and cr~(aeo+/9ei) = /3eo- In our model of the interaction 
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of a qubit with a given plant, we shall assume that the interaction is much 
shorter than r so that at most one qubit may interacting with a given plant 
at any instant of time. For two plants in cascade, we shall take them to be 
separated so that the time of flight of the qubits is exactly r seconds. This is 
purely for convenience and can be easily relaxed. For definiteness, we assume 
that each qubit is prepared independently in the "ground state" e and we 
denote by af the raising/lowering operators for the kth qubit: the operators 
corresponding to different qubits commute, while we have cr^T &t + a t a k = 1> 
(cr^) 2 = = (cr,r ) 2 . At time t\. = kr (k G N), we take the most recent qubit 
to interact with the first system to be the kth qubit, and the most recent to 
interact with the second to be the (k — l)st qubit. 

Let us denote the value of x > rounded down to the nearest whole 
number by [x\ and set 

of (*) := 

where a, (3 may take the values zero and one and where [B] u = 1 , [B] 1 = B for 
any operator B. In the following, we shall denote by O (r n ) any expression 
which is norm-convergent to zero as r — > as fast as r n . The identity 
rcr" 1 (A;) af (k) = af (k) + O (r) will be important in what follows and will 
correspond to the discrete version of the second order Ito products. For t > 
fixed, the processes 

Af (t) := r of (k) 
k=i 

are well-known approximations to the fundamental processes A a @ (t) in the 
limit r -> 0+, [9]. 

We shall fix bounded operators Hf on the jth system such that Hf^ = 

Hf and set H ( T j) (k) = Hf <8> af (k) . We shall first recall some well known 
results [9] for the situation where the qubits interact with only the first 
system (that is, set Hf = 0). The discrete time evolution is described 
by unitary kicks every r seconds according to U T (t) = U[t/ T \ ' " •U'M\ where 

Uk = exp |— WTi!^ (k)\. Expanding the exponential yields 14 = 1 + rGf (g) 

cr^ 13 (k) + O (r 2 ) with the Gf forming the coefficients of the unitary QSDE 

with parameters Gi related to Hi — \ H^l \ as in Appendix [A] 
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In the limit r — > + , the discrete time process U T (t) converges weakly in 
matrix elements to the solution of the QSDE 

dU (t) = Gf ® dA aP (t) U (t) . 

We now turn to the case of a cascaded system. This time the discrete time 
dynamics is given by V T (t) = V\t/ T \ • • • V 2 Vi where Vk = exp |— (k) — \tH^ (k — 1) 
Expanding the exponential now yields 

V k = 1 + rGf (g> of (k) + rGf (g> of (k - 1) + O (r 2 ) . 

with the Gf forming the coefficients of the unitary QSDE with parameters 
G2 related to H 2 as in Appendix [A] 

To better understand what is going on, we compute 

V fc V fc _x = 1 + rGf <g> of (k) 
+r jcf + Gf + GfG\ p \ ® of (k - 1) 
+rGf ® of (k-2) + (r 2 ) . 

This may be iterated to give 

V fc V fc _i--.V,= 

fc-1 

1 + r [Gf + Gf + Gf Gf } ® ^ < s (A; - 1) 

+rGf ® <xf (A;) + rGf ® cxf (/ - 1) + O (r 2 ) . 
Under the same mode of convergence as before, we obtain the limit QSDE 
dV t = G%- 1) ®dA a P (t) V(t) 

where we recognize Gf<_u = Gf + Gf + GfG^ 3 as the coefficients the uni- 
tary QSDE with the series product parameters G 2 < Gi, see fl39l) . Therefore 
G 2 ^i = G 2 < Gi. The generalization to multi-dimensional noise is straight- 
forward. 
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C Proof of Theorem 



5.6 



Clearly, if (1401) is satisfied, then both cascade systems are described by the 
same parameters, which implies that they are equivalent. Now suppose the 
two systems are parametrically equivalent, with S 2 undetermined. Now by 
Definition 15.31 we may obtain expressions for G2<Gi and Gi<G' 2 . Equating 
the first terms, we have S2S1 = S\S' 2 , and solving for S 2 one obtains S 2 = 
Sj^Si, as in (jUJ). Next, equating the second terms gives L 2 + S2L1 = 
Li + SiL' 2 . This expression can be solved for L/ 2 , as in (Hdl) . Similarly, the 
Hamiltonian term H' 2 in fill) can be found by equating the third terms. 

D Classical Systems as Commutative Quan- 
tum Subsystems 

In this subsection we explain how to model the classical system f lSTj) . shown 
in Figure [221, as a commutative subsystem of a larger quantum system. This 
representation is used in subsection 16.31 In equation floTI) . x{t) e R n , y{t) e 
R, f,g are smooth vector fields, h is a smooth real-valued function, and w 
and v are independent standard classical Wiener processes. 





classical system 









Figure 22: Block diagram of the classical system ( !B"T1 ). 



To model this classical system, we take the underlying Hilbert space of 
the system to be f) = L 2 (R n ) with q\ pj being the usual canonical position 
and momentum observables: q^ip (x) = Xjip (x) and pjijj (x) = —idjip{x). 
We write q = (q 1 , q n ) T , p = (pi, . . . ,p n ) T , and V = (<9i, . . . , d n ) T . If ip is 
a smooth function of x, then we find that, by Ito's rule, for tp t = (p(x(t)), 

dip = H c i ass i ca i (ip) dt + # T Vy? dw, (56) 

where £ c iassicai (v 9 ) — / T V</? + \g T ^ (g T ^ip) is the (classical) generator of the 
diffusion process x(t) in (J5TJ . 
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We seek a quantum network representation G c , as shown in Figure I2"3"l 













• 


At 


HD 






• 

L c -2 


A 2 





classical system 



Figure 23: Network representation of the classical system (1511) shown in 
Figure El 



The classical noises are viewed as real quadratures of quantum noises w = 
Ai + Al, v = A 2 + A 2 . Now define port operators L c \ = —ig T p— |V r g, L c2 = 
\h and internal Hamiltonian H c = \ (f T p + p T f) , where f = f — |[Vg]g 
(the Stratonovich drift) and g are n-vectors whose components are viewed as 
functions of q and h = h(q) is viewed as a self-adjoint observable function of q. 
We claim that the classical system (I5ip behaves as an invariant commutative 
subsystem of the open quantum system G c = (l,L cl ,if c ) ffl (l,L c2 ,0). To 
verify this assertion, we examine the dynamics. From (1311) we have 

dX c = (-i[X c , H c ] + C Lcl (X c ) + C Lc2 (X c ))dt 
+[X C , L cl ]{dA\ + dAi) + [X c , L c2 }(dA* 2 - dA 2 ) (57) 

Now set X c = (p = if(q), a, smooth function of the position operator. Then 
(|BT|) gives 

dtp = (-i[<p, H c ] + C Lcl (<p) + C Lc2 (<p))dt 

+[<p, L cl ] {dA\ + dA x ) + [if, L c2 ] (dA* 2 - dA 2 ) 

= (f T V V + l -g T V (g T V<p) )dt + g T V<p dw, (58) 

where, we have used -i[tp,H c ] = f T V(p, C Lcl (ip) = \g T V(g T Vtp), C Lc2 (ip) = 
0, [tp, L c i] = g T ^<£, and [<f, L c2 ] = 0. Hence the classical dynamics fl56|) is 
embedded in the dynamics of the position observable q only in the quantum 
system G q (independent of momentum dynamics). Note that only the real 
quadrature of the input field affects these dynamics, and they are unaffected 
by the field A 2 . 
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Next we look at the outputs. The first output is not of interest, so we 
focus on the second one. The output y{t) of the homodyne detector HD in 
Figure [23] is 

dy = dA 2 + dA* 2 = (L c2 + L* c2 )dt + dA 2 + dA\ = hdt + dv (59) 

which agrees with floTT) . as required. The unnormalized quantum filter for G c 
is 

da t (X c ) = a t (-i[X c ,H c ]+£ Lcl (X c )+jC Lc2 (X c ))dt 

+a t (L* c2 X c + X c L c2 )dy. (60) 

When X c = tp, this reduces to 

da t ((p) = cr< (£ c i ass i C ai (<p))dt + a t (h(p)dy, (61) 

which is the usual Duncan-Mortensen-Zakai equation of classical nonlinear 
filtering, [51 Chapter 18]. 
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